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ABSTRACT 

N 

The  singularity  expansion  method  (SEM)  has  been  applied  to 
determine  natural  resonances  of  a  horizontal  wire  and  perpendicular 
crossed  wires  oriented  over  an  Imperfect  ground  plane.  In  order  to 
account  for  the  imperfect  conductivity  of  the  ground,  the  Sotnmerfeld 
formulation  is  used  and  a  theoretical-numerical  solution  obtained. 
Sample  results  are  presented  for  both  the  frequency  domain  and  the 
SEM  solutions. 
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I.  INTRODUCTION 

An  electromagnetic  pulse  such  as  the  nuclear  EMP  or  the 
Lightning  EMP  will  induce  current  and  charge  densities  on  aircraft 
and  missile  surfaces.  These  distributions  represent  the  external 
electromagnetic  fields  that  are  in  general  related  linearly  to  the 
interior  fields  (equipment  currents  and  voltages).  Acquiring  a 
knowledge  of  the  exterior  fields  is  the  first  step  in  determining 
terminal  currents  and  voltages  Induced  within  an  illuminated  aero- 
nautical  system.  Transfer  functions,  either  theoretical  or  experi¬ 
mental,  are  used  to  relate  the  voltage  and/or  current  at  a  particular 
location  to  the  most  significant  point  or  points  of  entry  for  the 
electromagnetic  energy. 

A  number  of  simulation  facilities  have  been  constructed  to  test 
the  response  of  aircraft  and  missiles  to  the  nuclear  EMP.  In  each 
of  the  existing  facilities  the  aircraft  must  rest  on  a  concrete  pad 
(or  a  wooded  platform)  which  eliminates  the  possibility  of  testing 
the  Inflight  mode  response.  However  tests  can  be  conducted  as  the 

t 

aircraft  flies  by  the  simulator  but  this  procedure  presents  its 
own  difficulties.  Thus  a  theoretical  model  is  needed  for  the 
interpretation  of  the  "fly-by"  test  as  well  as  for  the  extrapolation 
of  the  test  data  from  ground  based  measurements  to  predict  the 
inflight  mode  response. 

In  contrast  to  nuclear-EMP  testing  llghtning-EMP  testing  has 
not  evolved  to  the  level  of  sophistication  of  nuclear  EMP  testing. 
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Probably  this  is  a  result  of  the  lightning  pulse  not  being  well 
defined  and  the  accompaning  nonlinear  environment  not  being  well 
understood.  However  llghtning-EMP  testing  requires  many  of  the 
same  considerations  as  nuclear-EMP,  simulating  the  appropriate 
electromagnetic  environment  and  quantifying  the  effects  of  the 
physical  limitations  of  the  simulator  configuration. 

In  order  to  gain  insight  into  the  response  of  an  aircraft 
or  missile  to  an  electromagnetic  pulse,  a  wire  model  is  used  for 
both  convenience  and  accuracy  fl}.  When  an  alrcraft/mlsslle  (or 
the  wire  model)  is  located  near  an  Imperfect  ground  plane  as  in 
the  usual  test  configuration,  the  induced  surface  currents  are 
affected  by  two  principle  processes.  First,  the  scatterer  is 
exposed  to  the  direct  radiation  and  the  ground  reflected  radiation. 
Second,  the  Induced  surface  currents  interact  with  the  ground  plane. 
The  first  process  is  well  understood  for  both  perfect  and  imperfect 
ground  planes,  and  is  straight  forward  to  analyze  [2].  However  the 
second  process  being  much  more  complicated  is  very  difficult  to 
analyze,  particularly  when  the  ground  is  an  Imperfect  conductor  [3]. 

A  general  formulation  based  on  the  singularity  expansion  method 
is  developed  for  horizontal  wire  scatterers  oriented  over  an  imper¬ 
fect  ground  plane.  In  order  to  account  for  the  Imperfect  conduc¬ 
tivity  the  exact  Somerfeld  formulation  is  used  [4].  Accordingly  a 
system  of  Integral  equations  are  derived  and  solved  utilizing  a 
numerical  solution  technique.  Natural  frequencies  for  a  single 
horizontal  wire  and  for  a  horizontal  wire  cross  are  obtained. 

Because  of  the  complexity  of  the  Somnerfeld  Integrals  as  they 


3 


Appeared  In  the  kernel  of  the  Integral  equations  the  computer  CPU 
tine  would  be  prohibitive  for  a  parametric  study.  Therefore 
sufficient  data  is  presented  only  for  fiducial  purposes.  Frequency 
domain  results  for  plane  wave  excitation  are  presented  along  with 
natural  frequencies. 
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II.  FORMULATION 

Frequency  Domain  Considerations 

For  an  arbitrary  configuration  of  horizontal  wires  oriented 
over  an  imperfect  ground  plane  the  individual  wire  currents  Induced 
by  an  Incident  electromagnetic  field  are  obtained  by  solving  a 
system  of  linear  Integral  equations  [5]. 


•  l  |  I1U'1)G(lj,£,i)dl’1  -  j4*«eo(Ej"C+  Er*f)  (1) 
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where  E*”C  and  E*®^  are  the  components  of  the  incident  and  ground 
reflected  electric  fields  along  the  jth  wire  of  a  system  of  N  wires, 
Il(*.'l)  is  the  current  on  the  ith  wire  at  position  I'j  , 
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Here  )  ,  the  coordinates  of  ^  along  the  ith  wire; 

(x  .y  .«,).  the  coordinates  of  points  ij  along  the  Jth  wire;  and 
the  angles  and'  are  Illustrated  in  Figure  1. 
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When  the  horizontal  wires  are  joined  additional  considerations 
must  be  made.  First  the  vires  must  satisfy  Kirchhoff's  current 
law,  the  sum  of  the  wire  currents  must  equal  zero  at  the  junction. 
Second  the  charge  per  unit  length  on  the  wires  are  required  to 
satisfy  the  Wu-King  [6]  junction  condition,  l.e.  for  an  n-wire 
junction 


’l*l  ■  ’2*2  '  ’3*3 . V. 

where  q^  is  the  charge  per  unit  length  given  by 

.  •  ’i  -  >  Iirti)l)„tta  pt 
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In  addition  to  the  junction  conditions  the  wire  currents  are 
required  to  vanish  at  the  open  ends. 

Due  to  the  complexity  of  the  system  of  integral  equations  to 
be  solved  a  numerical  solution  technique  is  employed,  in  particular, 
the  method  of  moments.  This  procedure  involves  dividing  the 
wire  structure  into  electrically  short  segments.  On  each  segment, 
the  wire  current  is  approximated  by  some  convenient  function,  here 
a  sinusoidal  current  expansion  is  used.  Therefore  the  current  on 
the  nth  segment  of  the  jth  wire  with  end  currents  I^Cl^.n)  and 
Ij(£j,n+1)  is  represented  by 


maUttiiiriiiaailaAdfcM 
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Upon  Introducing  the  segmental  current  representation  into  the 
system  of  Integral  equations,  a  system  of  linear  equations  for  the 
segment  end  currents  can  be  obtained  by  enforcing  the  resulting 
equations  at  a  discrete  set  of  points,  namely,  the  end  points  of 
the  wire  segments. 

The  resulting  system  of  linear  equations  is  of  the  form 

N 

l  s„i»  -  *.  <2» 

n*l 

where  N-N^+t^+N^ - is  the  total  number  of  unknown  currents 

after  requiring  the  currents  at  the  open  ends  to  be  zero.  The 

Integral  equation  is  enforced  at  the  ends  of  all  wire  segments 

excluding  open  wire  ends  and  junction  points,  which  yields  N-N 

w 

equations  where  is  the  total  number  of  wires  intersecting. 

Applying  the  Wu-King  junction  condition  yields  N  -NT  equations 
where  is  the  tota^  number  of  junctions.  Applying  the  Kirchhoff 

current  law  at  the  junctions  yields  an  additional  Nj  equations. 
Hence  the  total  number  of  equations  is  equal  to  N  the  total  number 

f 

of  unknown  currents. 

Computation  of  the  system  matrix  elements  S  can  be  simpll- 

ton 

fled  by  mathematical  manipulation  and  integration  by  parts.  How¬ 
ever  the  Sommerfeld  terms  accoutlng  for  the  Imperfect  ground  require 
a  double  numerical  integration,  one  Integral  over  the  wire  segment 


and  che  other  integral  over  an  infinite  range.  These  integrals 


occur  in  the  expressions  for  gcu  and  g  ,  see  (5)  and  (6).  To 

bit  V 

perform  the  evaluation  a  deformed  contour  similar  to  the  one 
suggested  by  Miller  et  al  [7]  is  chosen.  The  deformation  of  the 
contour  is  permlssable  since  all  the  branch  points  and  poles  of 
the  Integrand  lie  in  the  second  and  forth  quadrants  of  the  complex 
X  plane.  Accordingly 
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where  8*X+jX_  and  I  is  the  modified  Bessel  function  of  the 
X  o 

first  kind.  The  choice  of  X^  ,  a  real  constant,  as  the  limit  of 
the  first  integrals  in  (22)  and  (23)  is  made  to  render  the  inte¬ 
grands  of  the  second  Integrals  sufficiently  smooth  for  rapid 
convergence  when  a  Gaussian-Laguerre  interpolator^  quadrature 


formula  is  used.  For  real  frequencies  X^  •  5kQ/2ir  provides  good 
results.  Some  experimentation  is  required  to  choose  an  appropriate 
value  for  X^  when  complex  frequencies  are  conaidered.  The  first 
integrals  in  (22)  and  (23)  are  evaluated  by  using  the  Gaussian 
lnterpolatory  quadrature. 

Both  the  Gaussian  and  Gauss ian-Laguer re  quadrature  formulas 
are  derived  by  using  convenient  Interpolating  polynomials.  Theoret¬ 
ically,  the  more  roots  of  the  polynomials  (increasing  the  order) 
the  more  accurate  the  results  will  be.  However  this  occurs  at  the 
expense  of  computer  CPU  time.  For  the  data  that  will  be  presented 
subsequently,  the  number  of  roots  of  the  Gaussian  quadrature  is 
denoted  NIG  and  the  number  of  roots  for  the  GauBslan-Laguerre 
quadrature  is  NIL.  Both  NIG  and  NIL  are  varied  to  achieve  the 
desired  accuracy  within  the  limitations  of  available  CPU  time. 

Singularity  Expansion  Method 

In  order  to  employ  the  singularity  expansion  method  (SEM)  the 
foregoing  frequency  domain  formulation  is  extended  into  the  complex 
s-plane,  where  s  ■  j<o  .  Basically  the  SEM  solution  technique 
provides  a  solution  fop  the  Induced  current  in  terms  of  a  simple 
pole  expansion  in  the  frequency  domain  and  corresponding  damped 
sinusoids  in  the  time  domain.  To  construct  the  solutions  one  must 
obtain  the  natural  frequencies  (poles),  natural  modes  and  coupling 
coefficients  [8].  Only  the  natural  frequencies  are  expected  to 
sensitive  to  the  properties  of  an  imperfect  ground  plane.  Bence 
this  report  will  concentrate  on  obtaining  data  for  the  natural 
frequencies. 


The  natural  frequencies  are  the  simple  poles  of  the  solutions 
for  the  surface  current  and  charge  in  the  complex  s-domain.  In 
order  to  obtain  the  natural  frequencies  the  system  of  equations  in 
(21)  are  used.  By  observing  that  the  singularities  of  the  current 
are  those  complex  frequencies  for  which  the  system  matrix  is  slngu~ 
lar,  l.e.  det[S _ ]  ■  0  ,  the  natural  frequencies  are  obtained  by 

BUI 

searching  for  the  roots  of  the  determinant  of  the  system  matrix. 

For  the  data  reported  here  the  roots  were  obtained  via  a  Muller 
iteration  scheme  [9]. 
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III.  FREQUENCY  DOMAIN  RESULTS 


In  order  to  verify  the  numerical  procedure*  and  algorlthma 
frequency  domain,  domain  data  Is  obtained  for  twc  simple  wire 
configurations,  a  single  horizontal  wire  and  a  horizontal  wire 
cross.  Both  configurations  have  been  studied  extensively  by 
other  authors  for  both  free  space  conditions  and  perfect  ground 
plane  conditions  [10,11]. 

The  single  horizontal  wire  oriented  over  an  imperfect  ground 

a 

is  illustrated  in  Figure  2.  For  plane  wave  incidence  normal  to 
the  ground  and  with  the  electric  field  directed  parallel  to  the 
cylinder,  typical  induced  axial  currents  are  exhibited  in  Figure  3 
for  resonant  conditions  (for  comparison  see  Table  5).  Corresponding 
results  for  a  perfect  ground  would  show  significant  variation  of 
the  current  magnitude  with  height  above  the  ground.  From  the 
analysis  of  Taylor  et  al.  [12]  one  obtains  for  the  current  at  the 
center  of  the  horizontal  wire 


where 


T  J  4  U£lnc  r cos(ki/2)-l  1 

1  3  kZ  I  cos  (kit/2)  j 
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2n  2n 

u  “  n  +n  +  3k0*n  n  +n 

go  go 


(24) 


(25) 


(26) 


Figure  2:  Single  horizontal  cylinder  oriented 
over  an  imperfect  ground  plane. 
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to  the  cylinder. 


Ground  plane  ,  z  «  0 


A.  The  orthogonal  crossed  wire  configuration 
over  a  ground  plane. 


50  A/V/m 


Figure  5:  Current  distribution  induced  on 
perpendicular  crossed  cylinders 
oriented  horizontally  over  an 
imperfect  ground  plane  for  plane 
wave  illumination  propagating  normal 
to  the  ground  with  the  electric  field 
parallel  to  the  y-axis. 


IV.  SINGULARITY  EXPANSION  METHOD  RESULTS 


In  order  Co  have  confidence  in  the  numerically  obtained 
results  the  dependence  of  the  results  upon  the  various  numerical 
solution  parameters  is  investigated.  Considering  the  single 
cylinder.  Table  1  illustrates  typical  dependence  on  the  parameters 
used  to  evaluate  the  Sommerfeld  integrals  (22)  and  (23).  At  the 
lower  conductivity  substantial  variation  in  the  first  natural 
frequency  is  exhibited.  Also  the  natural  frequencies  depend  upon 
the  number  of  current  segments  used.  This  dependence  is  exhibited 
in  Table  2.  Generally  when  the  length  of  a  current  segment  is 
greater  than  a/z  and  less  than  2a,  good  results  are  obtained. 
Finally  the  numerically  obtained  results  are  compared  with  the 
results  of  independent  (but  less  rigorous)  formulations.  These 
data  are  shown  in  Table  3. 

Obtaining  the  natural  frequencies  for  various  ground  conduc¬ 
tivities  becomes  increasingly  more  difficult  as  the  conductivity 
is  decreased.  In  Figure  6  the  natural  frequencies  of  a  cylinder 

t 

are  displayed  in  the  complex  s-plane  with  conductivity  of  the 
ground  as  a  parameter.  For  high  conductivities  the  natural 
frequencies  approach  the  value  obtained  for  a  perfect  ground  [10). 
But  for  low  conductivities  numerical  difficulties  begin  to  occur 
for  0^  <  0.1  S/m  ,  which  prevented  convergence  to  the  free  space 

result.  The  dotted  line  is  expected  to  give  the  true  verletion. 
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TABLE  1:  First  Natural  Frequency  of  a  Cylinder 
over  an  Inperfect  Ground  for  Various 
Numerical  Parameters 


(t/a-20.,  h/a-4,  l- In,  c^-l) 


S1  i/c 


1.2S/m 


0.06S/m 


5k  /  2m 
o 

5k  I2it 
o 

5k  /2n 
o 

7k  /2n 
o 

9k  /2ir 
o 

Ilk  /2» 
o 


-0.1598  ♦  J  2.551 
-0.1598  ♦  J  2.551 
-0.1566  +  J  2.565 
-0.460  +J  2.584 
-0.442  +J  2.650 
-0.423  +  j  2.719 


TABLE  2:  First  Natural  Frequency  of  a  Cylinder 
over  an  Imperfect  Ground  Versus  the 
Number  of  Current  Segments 

(i/a-200,  h/a-40,  NIL -NIG -2, 

i-lm,  o-1.2xl08s/m,e  -1) 

R 


S1i/c 

-0.0293  + j  3.061 
-0.0392  +  «  3.018 
-0.0481  +J  2.980 
-0.05* *  +j  2.925 
-0.0570  +J  2.906 
-0.0593  +  J  2.878 
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TABLE  3:  First  Natural  Frequency  of  a  Cylinder 
over  an  Inperfect  Ground  as  Obtained 
by  Different  Analyses. 

(t/a-20,  h/a-4,  t-lm,  0-120  S/n  ,  eR-l  ) 


ANALYSIS 


Sommer f eld- Integral  Formulation 
Riggs  and  Shumpert  [13] 


S1£/c 

-0.0915  +j  2.599 
-0.0875  + j  2.526 


Reflection  Coefficient  Formulation  [14}  -0.1246  +  j  2.562 


TABLE  4:  First  Natural  Frequency  of  a  Cylinder 
over  an  Imperfect  Ground  versus 
Dielectric  Constant. 

(i/a-20. ,  h/a-4,  1-1. 0m.  0-0.12 S/m) 


S^/c 

-0.3320  +J  2.5182 
-0.3135  +J  2.5856 
-0.2741  +3  2.6140 
-0.2594  +J  2.6189 
-0.2379  +3  2.6232 


Note  that  when  numerical  difficulties  occur  the  natural  frequency 
Is  close  to  the  free  space  value. 

The  dependence  of  the  natural  frequencies  of  a  horizontal 
cylinder  upon  the  permittivity  of  the  ground  and  the  height  above 
the  ground  is  exhibited  in  Tables  4  and  5.  These  data  are  not 
intended  to  be  sufficient  to  provide  the  complete  parametric 
behavior.  They  are  presented  to  be  a  basis  of  comparison  for  a 
more  efficient  (l.e.  less  expensive)  algorithm  developed  by  using 
an  approximate  solution  technique. 
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TABLE  5:  First  Natural  Frequency  of  a  Cylinder 
over  an  Imperfect  Ground  Versus  Height. 

(l/a-20. ,  2.-1. Om,  o-0.12S/m,  e^-1) 

h/a  Sj^i/c 

4  -0.3320  4-  J  2.5182 

6  -0.2940  +  j  2.4899 

8  -0.3133  +  j  2.4324 

10  -0.3511  +  j  2.3768 

-0.4242  +  j  2.503 


TABLE  6:  Natural  Frequencies  of  Horizontal 
Crossed  Hires  over  an  Imperfect 
Ground  (i^tJ-L^l^  L/a-20,  ij/ij 

“0.5,  h/L-0.2,  L-lra,  o-0.1S/mt 
e  -20e  ) 

8  o 


Mode 

Free  Space 

SL/c 

Perfect  Ground 

Imperfect  Ground 

Sy  1,1 

-0.2923  + j 

2.319 

-0.0513  + j 

2.361 

-0.0965  + J 

2.222 

Sy  1,2 

-0.3426  + j 

3.726 

-0.1021  + J 

3.769 

-0.0927  +  J 

3. 715 

Sy  2,2 

-0.6786  + J 

6.06J6 

-0.3909  + J 

5.741 

-0.44  +j 

5.781 

sy  3,1 

-1.0166  +  j 

8.139 

-0.6691  + j 

7.581 

-0.69  +j 

7.59 

as  1 

-0.4242 +  j 

2.503 

-0.0898  + J 

2.592 

-0.2621  + J 

2.620 

V.  CONCLUSION 


By  extending  the  Sommerfeld  formulation  treating  Infinitesimal 
dipoles  over  an  imperfect  ground  plane,  a  formulation  is  developed 
for  treating  horizontal  wire  configurations  in  proximity  to  an 
Imperfect  ground.  Numerical  results  are  obtained  for  a  single 
wire  and  perpendicularly  crossed  wires. 

Structure  resonances  are  obtained  by  utilizing  the  singularity 
expansion  method.  Sample  results  are  presented  and  comparisons  are 
made  with  the  results  of  other  more  approximate  formulations. 
Because  of  the  lnordinant  amount  of  CPU  time  required  parameter 
studies  are  not  presented.  However  sufficient  data  are  provided 
to  verify  more  efficient  algorithms  that  may  be  developed. 
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